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We theoretically examine the effects of polaron formation in quantum dots on the 
transport properties. When a separation between two electron-levels in a quantum dot 
matches the energy of the longitudinal optical (LO) phonons, the polarons are strongly 
formed. The Rabi splitting between the levels is observable in a peak structure of the 
differential conductance G as a function of the bias voltage. The polaron formation 
suppresses the peak height of G, which is due to the competition between the resonant 
tunneling (resonance between a level in the dot and states in the leads) and the polaron 
formation (Rabi oscillation between two levels in the dot). G shows a sharp dip at 
the midpoint between the split peaks. This is attributable to the destructive interfer- 
ence between bonding and anti-bonding states in a composite system of electrons and 
phonons. 
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1. Introduction 

In semiconductor quantum dots, the relaxation of electrons between discrete energy-levels is 
an important issue from a viewpoint of application to quantum-dot lasers.^) The electron-phonon 
interaction should be much less efficient for the relaxation than in the bulk. This is because the 
number of phonons is very limited whose energies match the separation between the discrete lev- 
els,^' ^-^ which is in contrast to the case of bulk where the electron levels form a continuous spectrum. 
Related to this "phonon bottleneck," which is not always observed in actual quantum dots, several 
theoretical works have been reported on the relaxation mechanisms in quantum dots.^"^^ 

A remarkable exception of the inefficiency of the electron-phonon interaction is seen when the 
energy of the longitudinal optical (LO) phonons matches the level separation. Then the interaction 
can be significantly strong, reflecting a large density of states of the LO phonons with little disper- 
sion. The electron-LO-phonon interaction results in the formation of "polarons." The polarons in 
quantum dots are coherent states consisting of the discrete electron- levels and LO phonons, which 
were predicted by Inoshita and Sakaki.^'^^ The polarons have been found in self-assembled InAs 
quantum dots embedded in GaAs, using the far-infrared spectroscopy.^^' When a spacing be- 
tween the electron levels is tuned to match the energy of LO phonons by applying magnetic fields, 
the Rabi splitting of the coherent states has been observed in the absorption spectra. 

In the present paper, we theoretically examine a situation to observe the polarons by the 
electric transport, considering a quantum dot connected to external leads. The conductance takes 
place by the resonant tunneling through electron levels in the quantum dot. We show that the 
polaron formation is observable in a peak structure of the differential conductance G as a function 
of the bias voltage. Recently the transport properties of one of the self-assembled InAs quantum 
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dots can be investigated in vertically fabricated samples of heterostructures including a layer of 
such quantum dots.^^'^^^ Hence our theoretical results are expected to be observed experimentally. 
Compared with the optical experiments, the transport measurement has the advantage to be able 
to examine the polaron formation in a single quantum dot. Hence the measurement could reveal 
its characters in more details. 

Prom a viewpoint of the fundamental research, this is a simple example of important problems 
for the coupling between electrons and bosonic fields in the environment. We find new interference 
effects on the conductance in a composite system of electrons and phonons. (i) The polaron 
formation suppresses the peak height of G, which is due to the competition between the resonant 
tunneling (resonance between a level in the dot and states in the leads) and the polaron formation 
(Rabi oscillation between two levels in the dot), (ii) G shows a sharp dip at the midpoint between 
the split peaks. This is attributable to the destructive interference between bonding and anti- 
bonding states in the composite system. This is analogous to the transport through double quantum 
dots connected in parallel^^^ and Fano resonance. 

The organization of the present paper is as follows. In section 2, we explain our models and 
calculation methods. In section 3, we examine an exactly solvable model. By the calculations 
using Green's function, we show that the polaron formation is observable in a peak structure of 
the differential conductance. In section 4, we consider a more general model. The Green's function 
is calculated in the self-consistent Born approximation to take into account the electron-phonon 
interaction to infinite orders. The new interference effects on the conductance are discussed. The 
last section (section 5) is devoted to the conclusions and discussion. 

2. Models and Calculation Methods 

We consider a quantum dot with two levels, g and e, for electrons, as shown in Fig. 1(a). The 
level spacing, Ae = Se — Sg, is tuned. The levels are coupled to two external leads, L, R, through 
tunnel barriers. The Hamiltonian of electrons, Hg + i^x, is given by 

Hg = ^ei4di+ (1) 

i=g,e a=L,R k 

Ht = H ^(V^^4«a,fc + ll.C.), (2) 

i=g,ea=L,R k 

where c/|, aj^ ^ {di, a^^k) are the creation (annihilation) operators for an electron in the dot and 
in lead a, respectively. The spins of electrons are disregarded since they are not relevant in this 
problem. The electron-electron interaction is not taken into account in this paper. The tunnel 
couplings, i^T, broaden the level i in the dot with A' = -|- A^, 

al = hk:i', (3) 

where v is the density of states in the leads. 

For phonons, we consider an LO phonon mode, 

iJph = Y ^q^bq, (4) 

where 6q (bq) creates (annihilates) an LO phonon with momentum q. Since only phonons with long 
wavelength, |q|<27r/(dot size), interact with electrons significantly, we can regard the phonon 
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mode as dispersionless; Uq = cjlo- Other modes of phonons are neglected since they do not 
effectively couple to electrons in the quantum dot when Ae is of the order of lOmeV^^' (Appendix 
A). 

The elcctron-LO-phonon interaction is described by the Frohlich Hamiltonian i^c-ph (Appendix 
A). For the phonon-absorption process, the matrix element between an electron at level e with iV — 1 
phonon, |e, (A'" — l)q), and an electron at g with N phonon, \g, Nq), is written as 

(e, {N - l)q|i?e-ph|5, Nq) = ^qVq. (5) 

We denote 

which characterizes the strength of the electron-phonon interaction. We set v = 0.2^lO) consid- 
ering the experimental situations. 
The total Hamiltonian is 



Htot = He + Ht + -ffph + He-ph- (7) 

To examine the polaron formation, we have to include -ffe-ph to all orders. For the purpose, (i) 
we study an exactly solvable model in which only level e in the dot couples to the leads by the 
tunneling; = = 0. We present the Green's function in an analytical form in section 3. (ii) In 
general situations, we take Hq + + Hp\i as an unperturbed Hamiltonian, taking account of the 
tunnel couplings between the dot and leads exactly. Then we consider -ffe-ph as a perturbation, 
using the self-consistent Born approximation^^) (section 4). 



3. Calculated Results with Solvable Model 

To illustrate the effects of polaron formation on the transport properties, we begin with a 
simple model with V[ = = 0. This model is exactly solvable if we restrict the states in the 
quantum dot to |e,Oph) (an electron at level e and no phonon), 1^, Iq) (an electron at level g and 

one phonon), and no electron nor phonon (Fig. 1(b)). 
The Green's function G{£), which is defined by 

{£-Htot+i5)G{£) = l, (8) 

is obtained in an analytical way, as shown in Appendix B. The diagonal elements for the states in 
the dot are 

<^{e,Oph)(£) = (e,0ph|G(£)|e,0ph) 

= I (9) 
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and 



Gig,i,){e) = {g,lq\G{e)\g,l^) 



e-(£p + ^Lo)-S(g,i^)(e)' 


2 



In terms of the Green's function, the density of states in the dot is written as 



D{e) = --Im 

TT 



The electric current from lead L to i? is 



(11) 
(12) 

(13) 



2tt 



k,k' 



x[/(£fc-Mi) - fi^k' -/^r)], 



where fia is the Fermi level in lead a and /(s) is the Fermi distribution function. At T 
differential conductance is given by 

dl 



(14) 
0, the 



G = 



djiL 

-.1 



h {S-Se- ReE(e,Op,)(£))' + - Im^(e,Oph)(^)) 



(15) 



when /i/j is fixed. 

In Fig. 2, we show the calculated results of (a) the density of states in the quantum dot and (b) 
differential conductance G at T = 0. The level spacing, Ae = Se — eg, is changed gradually. In (a), 
two Lorentzian peaks appear corresponding to the states \g, Iq) and |e, Oph). The peak of the former 
{e ^ Eg + fkoi^o) is sharp and high, reflecting large density of states of the LO phonons. The peak 
of the latter (e ~ + Ae) is more broadened owing to the tunnel coupling to the leads. Around 
Ae = huJi^O: an anti-crossing between the states is seen, which is ascribed to the formation of 
polarons. The electron-phonon interaction mixes the states \g, Iq) and |e, Oph), making the bonding 
state with lower energy and the anti-bonding state with higher energy. Hence this anti-crossing of 
the two peaks is due to the Rabi splitting caused by the electron-phonon interaction. 

The differential conductance G shows a similar peak structure, as a function of the bias voltage, 
{Hl — fJ-R)/e, to that in the density of states D{e). The current flows by the resonant tunneling 
through the electronic states in the dot. Because of the mixture between \g, Iq) and |e, Oph), G has 
a peak even at the position of states \g, Iq) in spite of the absence of the tunnel coupling of level 
g. In conclusion, we clearly observe the Rabi-splitting by the polaron formation in the transport 
properties of the quantum dots. 



4. Calculated Results in Self-Consistent Born Approximation 



When both the levels in the quantum dot couple to the leads, we adopt the self-consistent Born 
approximation to consider the electron-phonon interaction to infinite orders (Appendix C). In this 
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approximation/^) the self-energy ^^(e) in Green's function of electrons, 



Gi(e) = \-. — — ^ (16) 



{i = g,e), is determined by the self-consistent equation 



Me) = ^ 



+ 



£ + hujq — Ej + iA^ — Tij (e -I- hcjq) 



(17) 



e — hujq — £j + lA^ — Sj(£ — hiOq)j 

Here, j ^ i and Nq is the number of phonons with momentum q}^^ Note that we disregard the 
inelastic transport processes where an electron enters a level in the dot and goes out from the other 
level in the dot, to focus on the effects of the polaron formation on the transport. This is a good 
approximation when ^ <C ^. Wc will discuss this approximation later. The density of states 
in the dot, D{e), and the differential conductance, G, are calculated in the same way as in the 
previous section. 

Figure 3 presents the calculated results at T = 0. The density of states in the dot shows three 
large peaks at e « Sg, £g + Hui^o, Sg + Ae, which correspond to 'level g with no phonon,' 'level g 
with one phonon' and 'level e with no phonon,' respectively. A small peak aX e ^ £g + Ae + 
is also observable, corresponding to 'level e with one phonon.' The peak widths arc determined 
mainly by the level broadening, A^'^, due to the tunnel couplings to the leads. An anti-crossing 
between \g, Iq) and le.Oph) is seen at Ae ~ ^lo? <is discussed in Fig. 2. 

The differential conductance G through the dot shows a similar peak structure as a function of 
the bias voltage. Now we discuss two new interference effects of the polaron formation on G. First, 
the polaron formation suppresses the peak height of G. When |A£ — ^lo| is large, the resonant 
tunneling through level g or e makes peaks of G ~ e^/h in height. When Ae ~ huji^o, however, 
G is much smaller than e'^/h. This is because the electron-phonon interaction results in a finite 
life-time of the electron levels, r (ImSj ~ —h/r), which weakens the resonant tunneling by a factor 
of ~ [A7(A^ + Vt)]^^°) 

Second, the differential conductance shows a dip at the midpoint between the two peaks. To 
see this dip more clearly, we show the conductance at T = with larger level broadening A^ in Fig. 
4. The level spacing in the dot is fixed to match the energy of the LO phonons; Ae = huj^o- This 
dip is sharp and deep when A*^ <C A*^. (In model (a) with A*^ = 0, the analytical solution indicates 
that G = exactly at jiL = £g + ^lo-) With increasing A^/A*^, the dip is smeared. This dip of 
G is a consequence of the destructive interference between the bonding and anti-bonding states of 
\g, Iq) and |e, Oph). Kawamura and Aono have proposed a similar dip of G in the transport through 
double quantum dots coupled in parallel.^^^ They have considered a situation where only one of 
the dots is connected to the leads by the tunneling. In this case, the differential conductance has 
two peaks due to the resonant tunneling through bonding and anti-bonding orbitals between the 
dots. G has a dip at the midpoint between the peaks (G = exactly when the tunnel couplings to 
left and right leads are equivalent). Another analogous phenomenon is the Fano resonance.^^^ The 
interference between a discrete level and a continuum of states results in an asymmetric shape of 
the resonance, with a dip on one side of the resonance. This dip is due to the destructive interaction 
between bonding and anti-bonding couplings of the discrete level and the continuum. Note that. 
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in our model, the dip of G is attributable to the destructive interference in a composite system of 
electrons and phonons. 

Finally, wc discuss the differential conductance G at finite temperature. When k-^T ^ A^'^, 
the peak heights of G are much smaller than /h, whereas the peak widths are determined by 
the thermal energy. Figure 5 shows G at k-QT = 0.05/kjlo {T ~ 18K), A£ = A|j = 0.005/kjlo 
and A| = A^ = O.OI/kjlo- Although the peaks of G are thermally broadened, the evidence of the 
polaron formation is still observable. We can also see the above-mentioned characters of G, the 
suppression of the peak heights at Ae ^lo and the dip of G between the resonant peaks. 

5. Conclusions and Discussion 

We have investigated the effects of polaron formation in a quantum dot on the transport 
properties. When a separation between two levels in the dot, g and e, matches the energy of the 
LO phonons, the polarons are strongly formed. The Rabi splitting between states \g, Iq) and |e, Oph) 
is observable in a peak structure of the differential conductance G, as a function of the bias voltage. 
We have found new interference effects on the transport in this composite system of electrons and 
phonons. The polaron formation suppresses the peak height of G, which is attributable to the 
competition between the resonant tunneling (resonance between a level in the dot and states in 
the leads) and polaron formation (Rabi oscillation between two levels in the dot). G shows a sharp 
dip at the midpoint between the split peaks. This is due to the destructive interference between 
bonding and anti-bonding states of l^, Iq) and |e, Oph). These theoretical results are expected to 
be observed experimentally. The measurement of the electric current would enable to observe the 
polaron formation in a single quantum dot and reveal its characters in more details. 

In the calculations with the self-consistent Born approximation, we have disregarded inelastic 
transport processes: When an electron is injected to a level in the dot, say e, from lead L, the 
electron is always ejected from the same level, e, to lead R. This is a good approximation when 
A£ ^ ^ A|^ ^. (Indeed the upper level usually couples to the leads more strongly than the lower 
level in quantum dots.) Otherwise, the relaxation from level e to ^ could take place in transport 
processes, which would broaden the peaks of the differential conductance. To consider the inelastic 
processes, we have to develop the self-consistent Born approximation in the nonequilibrium Green's 
function method. This is a challenging problem and requires the further study. Besides, we have not 
considered the coupling of electrons to longitudinal acoustic (LA) phonons^'^-' nor electron-electron 
interaction (Auger processes for the relaxation of electrons,^) etc.) which might cause the decay 
of polarons and, as a result, broaden the conductance peaks further. Although the quantitative 
estimation of these effects is beyond the scope of the present study, the polaron formation should 
not be influenced by them qualitatively. 

This work was partially supported by a Grant-in- Aid for Scientific Research in Priority Areas 
"Semiconductor Nanospintronics" (No. 14076216) of The Ministry of Education, Culture, Sports, 
Science and Technology, Japan. 
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Appendix A: Electron-Phonon Interaction in Quantum Dots 

The interaction between electrons and LO phonons is described by the Frohlich Hamiltonian.^^^ 
In quantum dots, it is given by 

^e-ph = -^J2M^'ij4dj{b.i + blg), (A-1) 

^ / ^ \ 1/4 

M„« = iyi5^-»(_) (i|e<---b), (A.2) 



/ m* / 1 1 



where \i) is the envelope function for electron level i in the dot and m* is the effective mass of 
electrons. £(oo) and e(0) are the dielectric constants at high and low frequencies, respectively. We 
assume that phonons are the same as those in the bulk GaAs {V is the volume of the bulk system) 
in which the InAs quantum dots are embedded. The dimensionless coupling constant cce-ph should 
be 0.15 to explain the experimental results, ^^'^^^ whereas a^-ph = 0.06 in bulk InAs. 

Owing to the factor of (i|e*''"^|j), the clcctron-phonon interaction is negligibly small when 
|q| > 2-11 /L with L being the dot size (~ 20nm^^'^^)). Hence we can regard the LO phonons as 
dispersionlcss; tiiOq = fkoi^o = 36meV. 

In this paper, we consider two electron levels (g,e). The coupling constant for the interlevel 
interaction is denoted by Vq = Mg^g^g/Vv in the text, Eq. (5). The strength of the couplings is 
characterized by 

2 \ ^ I i2 

V = \Vq\ 

|q|<27r/L 



|q|<27r/L 



2 



4«e-ph.j^ n2 / 1 /A 

= ^^(^LO)^^^^. (A.4) 

In the above-mentioned situation, i'/(^lo) ~ 0.2. We neglect the intralevel interaction, Mg^g^g, 
Mq^e,e) assuming that these effects are involved in the self-energies of electrons. 

We disregard the LA phonons for the following reasons, (i) LA phonons with \q\ > 2tt/L hardly 
couple to the electrons for the same reasons as LO phonons. (ii) LA phonons with \q\ < 27r/L are 
not relevant for the polaron formation since the energies (^2meV) are much smaller than the level 
spacing of electrons (Ae ~ lOmeV). 

Appendix B: Calculations of Green's Function (1) 

In the exactly solvable model with = = (Fig. 1(b)), we divide the total Hamiltonian 
as iftot = Hq + V, where Ho = He + i/ph and V = Ht + -ffe-ph- Using Eq. (8), 

(e, Oph|(£ -Ho + iS)G\e, Oph) = 1 + (e,Oph|FG|e,Oph), (B-1) 
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which yields 

{e-ee + i5)G(e,0ph) = 1 + E ^lid, UG\e, Oph) + Yl ^" (B-2) 

q a=L,R k 

Similarly, using Eq. (8), we obtain 

[e-{eg + fujq)+i6]{g,lq\G\e,0ph) = <G(e,Op,) (B-3) 
{e-ek + id){ak,0^\G\e,Opi,) = V^G^.^o^^y (B-4) 

Equations (B-2), (B-3) and (B-4) yield 



where 



ll/e|2 



a=L,H k 



-P ^ iT^5{£ - Ek) 

, e-Ek 



a=L,R 

-i{Al + A%), (B.6) 



and 



^(e,Oph)(^) - E 



^ E - {Eg + hWq) + iS 



t;2 



E- {Eg + fuJl,o) + iS 



(B-7) 



(Eqs. (9) and (10) in the text). The real part of St(£) can be included in the renormalization of 
the energy level £e, and hence it is omitted here. 

In the same way, wc obtain Xq)(£) analytical form of Eqs. (11) and (12). The electric 

current is written as Eq. (14) because 

(it;A:',Oph|r(£)|Lfc,Oph) = V'|*G(e,o,,)(£)F/, (B-8) 

where the T-matrix is given by T{£) = + H'yG{e)H'y. 

Appendix C: Calculations of Green's Function (2) 

When both the levels in the quantum dot couple to the leads, the Green's function cannot 
be obtained analytically. Then we divide the Hamiltonian as -fftot = ^^o + -f^c-ph with Hq = 
+ i^ph + i^x- First, we calculate the unperturbed Green's function, (e — Hq + iS)G^{£) = 1, in 
the same way as in Appendix B; 

G0(£) = (C-1) 

{i = g, e)}"^^ Next, the electron-phonon interaction, He-ph, is taken into account to infinite orders, 
by the self-consistent Born approximation.^^) The calculation yields Eqs. (16) and (17). 
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At r = 0, ATq = in Eq. (17). Sj(£) is expressed by Sj(£ — ^lo)) which is expressed in turn 
by Si(£ — 2hu!i,o), and so on. These yield the expression of ^^(e) as a continued fraction. We find 
that the continued fraction converges with the truncation at Ej(e — 6^lo)- At T 7^ 0, is 
coupled to T,i{e it hioi^o) in Eq. (12), Sj(e it hui^o) are coupled to T,i{e it 2?zu;lo) and Sj(e), and 
so on. These are simultaneous nonlinear equations. We solve the equations self-consistently in the 
approximation of Si(£ ± GHlolo) = to obtain Si(e). 

After the calculations of the Green's functions, we evaluate the density of states of electrons 
byi9) 

i.(£) = -ilm[G,(£)+Ge(£)]. (C.2) 
The electric current is expressed as 

^ = E E \Vh*Gi{ek)Vl\'6{ek' - £fc)[/(£fe - /xl) - f{ek' - tin)]- (C-3) 

1=3, e k,k' 

The differential conductance is obtained by G = edl/dfii. As mentioned in the text, we do not 
take account of inelastic transport processes where an electron enters a level in the dot and goes 
out from the other level in the dot. 
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(a) 




(b) 



V 



|e, Oph) 



Fig. 1. (a) We consider a quantum dot with two levels, g and e, in the presence of electron-LO-phonon interaction. 
The level spacing, Ae = Se — £g, is tuned. The levels are coupled to two external leads {a — L, R) through tunnel 
barriers with and V^. (b) In a solvable model with Vl — = 0, the states in the dot are restricted to |e,Oph) 
(an electron at level e and no phonon; middle panel), \g, Iq) (an electron at level g and one phonon; right panel), 
or no electron nor phonon (left panel). 



Tasai and Eto (Fig. 1) 
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Fig. 2. The calculated results of an exactly solvable model with Vl = — 0. The broadening of level e by the 
tunnel coupling to the leads is Al — A|j = O.OlfitJLO- (a) The density of states in the quantum dot, D{e), with 
changing the level spacing Ae = ee — Eg. (b) The differential conductance G through the dot at T = 0, as a function 
of the Fermi level in lead L, hl, whereas fiR is fixed. 



Tasai and Eto (Fig. 2) 
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(a) 



(b) 




Fig. 3. The calculated results using the sclf-consistont Born approximation. The broadenings of levels g and e by the 
tunnel coupling to the leads are = Af^ = O.OOS^icjlo and A1 — A% = O.Oltiuii^o, respectively, (a) The density 
of states in the quantum dot, D(e), with changing the level spacing Ae — e^ — £g- (b) The differential conductance 
G through the dot at T = 0, as a function of the Fermi level in lead L, hl, whereas hr is fixed. 



Tasai and Eto (Fig. 3) 
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Fig. 4. The differential conductance G through the quantum dot at T = 0, as a function of the Fermi level in lead L, 
Hl ifJ-R is fixed). The level spacing in the dot matches the energy of the LO phonons; Ae — huji^o- The broadening 
of the upper level e is fixed at A£ = = 0.2/t(j;LO, whereas that of the lower level g is = A|j = 0.02ftajLO 
(solid line) and 0.05/iu;LO (dotted line). The case without the tunnel coupling of the lower level (Vl = V^ =Q) is 
shown by broken line. 



Tasai and Eto (Fig. 4) 
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Fig. 5. The differential conductance G through the quantum dot at ksT — 0.05rt(j;LO {T ~ 18K), as a function of 
the Fermi level in lead L, jj,L {hr is fixed). The level spacing As — — £g in the dot is gradually changed. The 
level broadenings are the same as in Fig. 3 (A^ = = 0.005/ia;LO, A|, = A|j = O.OIRwlo)- 



Tasai and Eto (Fig. 5) 



